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Abstract 

The  phenomena  of  superheating  and  supercooling  can  lead  to  a  thermal  hysteresis  curve  where  two  states  are  stable 
over  a  range  of  temperatures.  Here  we  present  a  simple  magnetic  multilayer  system  (Fe/Gd)  that  can  be  designed  to 
show  a  thermal  hysteresis  curve.  Calculations  show  that,  with  proper  choice  of  parameters,  the  width,  in  temperature,  of 
the  hysteresis  curve  can  be  controlled  by  an  external  magnetic  field  and  varies  from  20  to  100K  over  a  field  range  of 
300-600  Oe.  Polarized  neutron  reflectivity  measurements  confirm  that  such  behavior  is  realized  experimentally  in 
multilayers  of  Fe/Gd. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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Superheating  and  supercooling  are  well-known 
phenomena.  Most  of  the  examples  of  superheating 
and  supercooling  concern  physical  changes  of  state 
such  as  melting  or  crystallization  [1].  A  magnetic 
example  of  supercooling  has  been  discussed 
recently  regarding  magnetic  phase  transitions  in 
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erbium  [2].  In  magnetism  one  can  also  find 
discussions  of  high-  and  low-spin  transitions  in 
molecules  that  have  an  associated  thermal  hyster¬ 
esis  [3],  Many  of  these  systems  have  a  requirement 
of  ultrapure  samples  since  impurities  can  act  as 
nucleation  sites  and  initiate  the  phase  transition. 

Here  we  present  a  simple  magnetic  multilayer 
system  that  can  be  designed  to  show  a  thermal 
hysteresis  curve.  The  width,  in  temperature,  of  this 
hysteresis  curve  can  be  controlled  by  an  external 
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magnetic  field  with  small  fields  (200-300  Oe) 
leading  to  thermal  hysteresis  spanning  about  a 
hundred  degrees  Kelvin.  Larger  fields  (600  Oe  and 
above)  reduce  the  width  of  the  hysteresis  curve  to 
20  K. 

We  consider  a  multilayer  system  with  alternat¬ 
ing  ferromagnetic  films  that  are  antiferromagneti- 
cally  coupled  at  the  interfaces.  In  addition  it  is 
necessary  that  the  magnetization  of  one  material 
change  more  rapidly  with  temperature  than  the 
other.  An  example  of  such  a  system  is  the  Fe/Gd 
multilayer  where  Gd  thermal  averaged  magnetic 
moment  changes  from  7  /tB  to  zero  as  T  changes 
from  0  to  about  300  K.  There  are  a  number  of 
known  magnetic  states  for  the  Fe/Gd  system  [4-6] 
which  are  illustrated  in  Fig.  1.  There  is  a  low- 
temperature  aligned-Gd  state  where  the  Gd 
magnetization  is  aligned  with  the  external  field 
(and  the  Fe  is  opposite)  and  there  is  a  high 
temperature  aligned-Fe  state  where  the  Fe  magne¬ 
tization  is  aligned  with  the  external  field.  These 
configurations  are  stable  when  the  applied  mag¬ 
netic  field  is  small;  in  higher  magnetic  fields  there  is 
a  spin-flop  like  state  where  the  Fe  and  Gd 
moments  are  both  canted  with  respect  to  the 
applied  field. 

The  physical  origin  for  the  thermal  hysteresis  is 
simple.  At  high  temperatures  the  system  is  in  the 
Fe-aligned  state;  the  Gd  moments,  which  are 
small,  are  oppositely  directed  to  the  field.  As  the 
temperature  decreases  the  Gd  moments  increase 
and  the  net  magnetic  moment  can  actually  be 
opposite  to  the  external  field  if  an  anisotropy  in  Fe 
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holds  the  Fe  spins  in  place.  Eventually,  however, 
this  configuration  becomes  unstable  and  the 
structure  reverses  with  the  Gd  moments  along 
the  field  and  the  Fe  moments  antiparallel  to  the 
applied  field.  A  similar  situation  occurs  on  heating 
the  system  from  low  temperatures. 

The  general  situation  of  a  Fe/Gd  multilayer  is 
quite  complicated  and,  in  fact,  shows  dramatic 
surface  effects,  as  we  shall  see.  We  can,  however, 
obtain  some  useful  intuition  using  a  simple  model 
which  includes  the  Zeeman  energy  of  each  film  in 
the  applied  field,  the  interfacial  exchange  energy 
between  the  two  films  and  a  uniaxial  anisotropy  in 
film  1  (Fe). 

E  =  —  H  *  (t\m\  cos  6 1  +  t2m2  cos  (L) 

—  Jm\m2  cos (02  —  0i)  —  ( 1  w  1  H.d cos2  0 (1) 


Here  H  is  the  external  magnetic  field,  t\  the 
thickness  of  film  1,  mx  the  magnetization  of  film  1 
and  0i  is  the  angle  between  ni\  and  the  external 
field.  Similar  definitions  hold  for  the  2nd  material. 
Ha  is  the  anisotropy  field,  assumed  to  exist  in  film 
1  (Fe)  only  for  simplicity,  and  assumed  to  be  in  the 
same  direction  of  the  measuring  field.  We  note  that 
Gd  is  in  an  S  state  and  therefore  the  anisotropy  is 
expected  to  be  small.  In  practice  anisotropy  can  be 
induced  by  different  types  of  growth  methods 
(deposition  in  a  magnetic  field,  tilting  the  substrate 
during  deposition,  etc.).  However,  the  general 
features  obtained  here  should  persist  if  the 
anisotropy  is  found  in  the  Fe,  in  both  materials, 
or  only  in  the  Gd.  J  measures  the  exchange 
coupling  between  films  and  is  negative  for 
antiferromagnetic  coupling. 

To  have  a  stable  magnetic  state  we  require  the 
usual  conditions  that  £(01,  02)  is  a  minimum 


0£  _  0£ 


and 
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Fig.  1.  Illustration  of  the  three  main  phases  in  the  Fe/Gd 
superlattice.  In  the  Gd-aligned  phase  the  thermal  averaged 
moments  of  the  Gd  are  large  and  aligned  with  the  external  field. 
In  the  Fe-aligned  phase  at  higher  temperature  the  Gd  moments 
are  small  and  the  Fe  moments  align  with  the  field. 
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It  is  easy  to  show  that  Eqs.  (2)  and  (3)  are  satisfied 
by  either  the  Fe-aligned  state  ($i  =  O,02  =  ft),  or  by 
the  Gd-aligned  state  (0l=n,d2  =  0),  as  long  as  the 
external  field  is  small  compared  to  the  exchange 
field.  The  critical  condition  of  Eq.  (4)  can  be  found 
using  Eq.  (1)  and  solved  for  m2,  the  moment  in  the 
Gd  film.  After  some  algebra,  we  find 

ti(±H2ti  +  mxHJ  +  2 H.AHh  ±  2HdmxJ) 

,U2  = - " - JhH - " - ' 

(5) 

Elere  the  +  sign  indicates  the  critical  value  of  m2  for 
the  high-temperature  state  (Fe-aligned)  and  the  — 
sign  indicates  the  critical  value  for  the  low- 
temperature  state  (Gd-aligned). 

We  show  the  stability  limits  from  Eq.  (5)  in  Fig. 
2.  The  parameters  for  the  calculation  are  chosen  to 
reasonably  match  the  experimental  conditions 
discussed  below.  The  thickness  of  the  Fe  film  is 
?i  =  35x  10_8cm;  the  thickness  of  the  Gd  film  is 
t2  =  50  x  10_8cm;  the  uniaxial  anisotropy  in  Fe  is 
Hd  =  50  G.  The  magnetization  in  Fe  is  assumed  to 


Fig.  2.  Stability  limits  for  the  Fe-  and  Gd-aligned  states  as  a 
function  of  applied  magnetic  field.  At  low  applied  fields  both 
states  are  stable  over  a  range  of  magnetization  values  for  Gd. 
(a)  shows  the  phase  diagram  for  a  larger  antiferromagnetic 
interfacial  coupling  and  (b)  shows  the  phase  diagram  for 
smaller  antiferromagnetic  interfacial  coupling.  As  the  antifer¬ 
romagnetic  coupling  is  increased  the  region  where  both  states 
are  stable  also  increases. 


be  constant  over  the  temperature  range  measured 
and  is  given  by  ni\  =  1740 G.  Thus  this  figure 
shows  the  values  of  the  Gd  magnetization  at  which 
the  different  states  become  unstable.  For  example, 
at  low  applied  field  the  Gd-aligned  state  is  stable 
for  m2  values  above  the  lower  curve,  and  the  Fe- 
aligned  state  is  stable  for  m2  values  below  the 
upper  curve.  There  is  also  a  range  of  Gd 
magnetization  values  where  both  states  are  stable. 
This  range  of  magnetization  values  corresponds  to 
a  range  of  temperatures  since  the  Gd  magnetiza¬ 
tion  changes  rapidly  with  temperature.  Because 
there  are  two  stable  states  in  this  temperature 
range,  one  has  the  possibility  of  thermal  hysteresis. 
As  the  applied  field  is  increased,  this  region  where 
the  two  states  are  simultaneously  stable  becomes 
much  smaller.  This  shows  that  the  thermal 
hysteresis  curve  ought  to  narrow  substantially  at 
higher  fields.  In  fact,  by  setting  the  two  critical 
values  in  Eq.  (5)  equal  to  each  other,  we  can  obtain 
an  equation  for  the  maximum  applied  field  which 
allows  thermal  hysteresis.  We  find 

Hm  =  y/-2  HamxJ/t2.  (6) 

A  comparison  of  Fig.  2(a)  and  (b)  shows  the 
effect  of  interfacial  exchange  on  the  stability 
regions  for  the  different  states.  The  most  notice¬ 
able  feature  is  that  when  the  exchange  coupling  is 
weaker,  the  region  where  both  aligned  states  are 
stable  is  restricted  to  a  smaller  range  of  applied 
field  values.  The  reason  for  this  is  that  a  weaker 
interfacial  exchange  allows  the  canting  of  spins  at 
the  interface  at  lower  applied  fields.  The  aligned 
states  become  unstable  and  the  spin  flop  state  is 
the  only  one  allowed.  Similarly,  a  reduction  of  the 
anisotropy  also  reduces  the  range  for  thermal 
hysteresis  in  agreement  with  Eq.  (6). 

We  can  make  a  direct  connection  between  the 
results  in  Fig.  2  and  the  temperature  range  over 
which  both  states  are  stable  by  using  a  graph  for 
magnetization  as  a  function  of  temperature  for  Gd 
as  shown  in  Fig.  3.  For  example,  at  H  —  200  Oe  in 
part  (b)  of  Fig.  2  we  see  that  both  states  are  stable 
when  the  Gd  magnetization  is  between  0.73  and 
1.47  kG.  Using  Fig.  3,  we  can  see  these  magnetiza¬ 
tion  values  correspond  to  temperatures  of  272  and 
200  K,  a  range  of  72  K. 
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The  calculation  associated  with  Eqs.  (1)  (5) 
assumes  that  the  spins  within  a  him  are  rigidly 
locked  together.  If  the  exchange  coupling  between 
films  is  significant,  this  may  not  be  a  good 
assumption.  We  therefore  examine  the  thermal 
hysteresis  in  the  Fe/Gd  system  using  a  self- 
consistent  local  mean-held  method  that  has  been 
used  successfully  by  a  number  of  groups  to 
describe  this  system  [4-7]. 

In  this  model  we  consider  alternate  hlms  of  BCC 
Fe  and  HCP  Gd.  The  spins  in  each  monolayer  are 
ferromagnetically  ordered  and  lie  parallel  to  the 
him  surfaces.  Only  nearest-neighbor  exchange 
interactions  are  considered.  A  spin  in  monolayer 
i  feels  a  total  held  composed  of  the  exchange  helds 
from  the  layers  above  and  below,  external  helds 
and  anisotropy  helds: 

h,  =  <£■+!>  +  Jij-1  (Si- 1» 

nn 

+  m+H,(S*)/Si.  (7) 

We  use  an  anisotropy  of  50  G  in  Fe  and  Ha  is 
taken  to  be  zero  in  Gd.  All  other  parameters  are 
given  in  Ref.  [4].  (The  reduction  parameter  of  Ref. 
[4]  is  0.75.).  The  minimum  energy  state  is  found 
using  an  iterative  method  where  a  layer  of  spins  is 
chosen  at  random.  One  hnds  the  net  effective  held 
acting  on  this  layer  and  the  spins  in  the  layer  are 
rotated  to  lie  along  this  held.  This  lowers  the 
energy  of  the  system.  The  process  is  repeated  until 
a  self-consistent  state  is  reached  where  the  spins  in 


17  layers  Gd  /  23  layers  Fe 


Temperature  (K) 

Fig.  4.  Theoretical  thermal  hysteresis  curves  for  a  17  ML  Gd/ 
23  ML  Fe  superlattice  for  different  applied  magnetic  fields.  The 
hysteresis  curves  initially  narrow  as  the  applied  field  is 
increased.  At  higher  fields  (12000e)  the  system  is  in  a  spin- 
flop  phase  and  does  not  show  substantial  hysteresis. 

each  layer  are  aligned  with  the  effective  helds 
produced  by  the  external  held,  the  anisotropy 
helds  and  the  exchange  helds  of  the  neighboring 
layers.  To  take  temperature  into  account,  the 
thermal  averaged  magnitude  <S,)  is  given  by  the 
Brillouin  function  <  -S',)  =  SiB(S,Hj/kT). 

The  results  for  the  Fe  magnetization  as  a 
function  of  temperature  for  an  infinite  Fe/Gd 
superlattice  are  given  in  Fig.  4.  This  curve  is 
calculated  for  a  structure  with  17  monolayers  Gd 
and  23  monolayers  of  Fe  in  each  unit  cell.  Again 
this  choice  of  values  is  intended  to  match  the 
experimental  parameters  of  Gd  50  A/Fe  35  A.  (We 
assume  c-axis  growth  for  the  Gd  with  a  lattice 
parameter  of  5.78  A  and  two  layers  per  unit  cell; 
for  Fe  the  layer  to  layer  distance  is  about  1.43  A). 
The  uniaxial  anisotropy  in  Fe  is  taken  as  50  G. 
Periodic  boundary  conditions  are  used  to  simulate 
the  inhnite  superlattice.  The  width  of  the  hysteresis 
curve  is  quite  large  for  low  magnetic  helds  and 
becomes  narrower  at  higher  helds  in  agreement 
with  the  results  of  Fig.  2.  For  even  higher  helds  the 
system  does  not  display  hysteresis,  but  makes  a 
transition  from  the  aligned  phases  to  the  spin  hop 
phase  and  the  magnetization  essentially  increases 
linearly  as  a  function  of  temperature.  Again  this  is 
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in  agreement  with  Fig.  2  where  we  have  seen  that 
neither  the  Gd-aligned  nor  the  Fe-aligned  state  is 
stable  in  the  high-held  region  of  the  phase 
diagram. 

It  is  known  that  the  finite  size  of  a  Fe/Gd 
superlattice  can  have  a  significant  influence  on  the 
phase  transitions.  For  example,  at  low  tempera¬ 
tures  and  low  fields  the  system  is  typically  in  the 
Gd-aligned  state.  If  the  magnetic  held  is  increased 
the  magnetic  structure  will  eventually  change  to  a 
spin-hop-like  state.  Flowever,  the  critical  held  at 
which  this  phase  transition  takes  place  depends 
dramatically  on  whether  the  outermost  surface  is 
Gd  terminated  or  Fe  terminated.  If  the  structure 
has  Fe  spins  on  the  outside,  then  in  the  Gd-aligned 
state  the  Fe  spins  are  opposite  to  the  external  held 
and  a  surface  phase  transition  to  the  spin  hop  state 
nucleates  at  the  Fe  surface  [8-10],  The  critical  held 
for  this  surface  transition  is  about  a  factor  of  5 
lower  than  that  for  the  bulk  transition. 

The  thermal  hysteresis  effect  discussed  here  also 
shows  substantial  hnite  size  effects.  In  Fig.  5,  we 
compare  results  for  an  inhnitely  extended  multi¬ 
layer  with  a  unit  cell  17  ML  Gd/23ML  Fe  with 
those  of  a  hnite  multilayer  with  the  same  unit  cell 
but  with  a  total  of  15  repetitions.  The  thermal 
hysteresis  for  an  inhnite  structure  is  approximately 
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Fig.  5.  Theoretical  thermal  hysteresis  curves  for  a  finite  and 
infinite  Fe/Gd  multilayer.  The  unit  cell  is  17  ML  Gd/23  ML  Fe 
in  both  cases.  The  applied  field  is  300  Oe.  We  note  that  the  finite 
multilayer  has  a  substantially  narrower  thermal  hysteresis 
curve. 


twice  as  wide  as  that  for  the  hnite  structure.  This  is 
due  to  the  surface  phase  transition  discussed 
above.  In  a  hnite  structure  composed  simply  of  a 
repetition  of  the  Fe/Gd  unit  cell,  the  spins  at  one 
surface  will  always  be  opposite  to  the  applied  held 
in  an  aligned  state,  and  as  a  result  a  surface  phase 
transition  will  nucleate  at  that  surface.  These 
surface  effects  will  be  discussed  in  detail  elsewhere. 
Nonetheless,  the  main  features  of  the  thermal 
hysteresis  remain  qualitatively  the  same  regardless 
of  the  number  of  unit  cells,  and  it  is  important  to 
compare  our  theoretical  results  with  the  results  of 
an  experiment.  In  fact,  in  a  related  system,  a 
thermal  hysteresis  curve  has  already  been  seen  in  a 
reorientation  transition  [1 1]  of  ultrathin  Fe  on  Gd. 

Measurements  were  taken  on  a  multilayer  with 
composition  [Gd  50A/Fe  35  A]  x  15  deposited  on 
a  Si  wafer  with  Fe  as  the  top  layer  of  the 
multilayered  stack.  The  sample  was  covered  with 
40  A  Si  to  prevent  oxidation.  X-ray  rehectometry 
conhrms  a  well  ordered  sample  with  hat  interfaces 
(rms-roughness  7  A).  The  compensation  tempera¬ 
ture  was  determined  to  be  125  K  by  Squid 
magnetometry.  The  magnetic  properties  were 
studied  with  polarized  neutron  rehectometry  at 
the  Intense  Pulsed  Neutron  Source  of  Argonne 
National  Laboratory.  Fig.  6  presents  the  intensity 
of  the  hrst  Bragg  reflection  of  the  superlattice  as  a 
function  of  the  neutron  momentum  transfer 


Fig.  6.  The  spin-dependent  neutron  reflectivities  in  the  region 
of  the  first  Bragg  reflection  for  the  two  polarization  states  of  the 
neutron.  R+  (black  circles)  is  the  reflectivity  for  the  case  where 
the  incident  neutrons  are  polarized  parallel  to  the  applied  field 
and  R  (open  circles)  refers  to  incident  electrons  polarized 
antiparallel  to  the  applied  field. 
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perpendicular  to  the  surface,  for  the  two  polariza¬ 
tion  states  of  the  neutron  [12].  The  width  of  the 
peak,  which  is  inversely  proportional  to  the 
number  of  layers  coherently  diffracting  neutrons, 
as  well  as  the  presence  of  secondary  oscillations, 
due  to  the  total  superlattice  thickness,  indicate  the 
high  perfection  of  the  sample. 

In  the  Fe  aligned  state  the  intensities  diffracted 
at  the  first  Bragg  peak  are  proportional  to  a 
combination  of  the  scattering  amplitudes  of  the 
single  bilayer: 

/±  OC  [( N Fe^Fe  ~  Aod^Gd)  ±  (JA^Fe^Fel 

+  I^Gd^Gdl)]2,  (8) 

whereas  for  Gd  aligned  in  field  direction  the 
intensities  are  proportional  to: 

E  OC  [(-IV Fe^Fe  —  -(Vod&Gd)  T  (I^Fe/’Fel 

+  I^Gd^Gdl)]2,  (9) 

Z?Fe,Gd  are  the  nuclear  and  pFc.c,d  the  magnetic 
scattering  lengths  for  the  two  elements  that  are 
assumed  aligned  in  the  direction  of  the  applied 
field.  A  convenient  quantity  is  the  spin  asymmetry 
P  —  (/+  —  I~)/(I+  + 1~),  which  is  basically  pro¬ 
portional  to  the  magnetization  of  the  Fe  sublattice, 
and  thus  it  can  be  immediately  compared  with  the 
theoretical  predictions.  Fig.  7  shows  the  spin 
asymmetry  when  the  sample  is  cycled  as  a  function 
of  temperature  at  different  fields.  For  H  <  1200  Oe 
the  loop  shows  hysteresis.  The  area  subtended  by  the 
loop  decreases  as  the  field  is  increased;  yet  the  slope 


Fig.  7.  Thermal  hysteresis  of  the  asymmetry  P.  The  vertical  axis 
is  proportional  to  the  Fe  magnetization. 


of  the  segments  uniting  the  lower  to  the  upper 
branch  remains  the  same.  Above  900  Oe  the  loop  has 
closed  down  almost  completely,  but  with  a  smoother 
transition  from  the  lower  to  the  upper  branch. 

The  results  presented  here  show  the  excellent 
qualitative  agreement  of  theory  and  experiments, 
particularly  in  that  the  width  of  the  thermal 
hysteresis  decreases  as  the  external  field  is  in¬ 
creased,  and  in  the  transition  to  the  spin  flop  state 
at  higher  fields,  a  state  where  the  magnetization 
increases  approximately  linearly  with  temperature. 

It  would  be  desirable  to  have  a  direct  experi¬ 
mental  confirmation  that  at  low  field  the  magne¬ 
tization  simply  reverses,  while  at  higher  fields  the 
magnetization  reverses  by  gradual  rotation,  from 
an  aligned  state  to  a  spin-flop  configuration  to  a 
reversed  aligned  state.  This  verification  was  not 
possible  in  our  experiments  due  to  technical 
reasons.  During  reversal,  the  system  breaks  down 
in  domains  of  ~  1000  A  lateral  dimension.  These 
cause  broadening  of  the  reflected  beam.  The  total 
intensity  is  recorded  accurately  by  the  detector, 
but  the  used  supermirror  polarization  analyzer 
operates  correctly  only  on  well-collimated  beams. 
An  experimental  check  with  a  3He-polaization 
analyzer  which  accepts  a  broad  angular  spread  of 
the  reflected  beam  would  be  desirable  to  check  the 
spin-flop  transition.  However,  even  without  that 
verification,  abundant  evidence  for  this  transition 
already  exists  [13-18]. 

In  summary  we  have  demonstrated  both  theo¬ 
retically  and  experimentally  a  novel  system  which 
has  a  thermal  hysteresis  curve  that  can  be 
significantly  changed  by  application  of  an  external 
magnetic  field.  Although  the  present  work  deals 
with  a  specific  system,  Fe/Gd,  it  is  applicable  to  a 
whole  set  of  layered  structures  where  the  tempera¬ 
ture  dependence  of  one  material  is  substantially 
different  from  that  of  the  other. 


Acknowledgments 

The  work  by  REC  was  supported  by  a  grant 
from  the  ARO  Grant  IDAAD 19-02- 1-0 174.  The 
work  by  GPF  was  supported  by  US  DOE,  Office 
of  Science  Contract  W31-109-ENG-38.  WL  ac- 


R.E.  Camley  et  al.  /  Journal  of  Magnetism  and  Magnetic  Materials  286  (2005)  65-71 


71 


knowledges  financial  support  by  NATO  CRG 
programme  SA  (CRG.CRG.973129). 


References 

[1]  Z.H.  Zhang,  P.  Kulatunga,  H.E.  Elsayed-Ali,  Phys.  Rev. 
Lett.  56  (1997)  4141. 

[2]  K.A.  Gschneidner  Jr.,  V.K.  Pecharsky,  David  Fort,  Phys. 
Rev.  Lett.  78  (1997)  4281. 

[3]  O.  Kahn,  C.  Jay  Martinez,  Science  279  (1998)  44. 

[4]  M.  Sajieddine,  Ph.  Bauer,  K.  Cherifi,  C.  Dufour,  G. 
Marchal,  R.E.  Camley,  Phys.  Rev.  B  49  (1994)  8815. 

[5]  See  the  review  article  by  R.E.  Camley,  R.L.  Stamps,  J. 
Phys.:  Condens.  Matter  5  (1993)  3727. 

[6]  R.E.  Camley,  D.R.  Tilley,  Phys.  Rev.  B  37  (1988)  3413; 
R.E.  Camley.  Phys.  Rev.  B  39  (1989)  12316. 

[7]  K.  Takanashi,  Y.  Kamiguchi,  H.  Fujimori,  M.  Motokawa, 
J.  Phys.  Soc.  Jpn.  61  (1992)  3721. 

[8]  J.G.  LePage,  R.E.  Camley,  Phys.  Rev.  Lett.  65  (1990) 
1152. 


[9]  D.  Haskel.  Y.  Choi,  D.R.  Lee,  J.C.  Lang,  G.  Srajer.  J.S. 
Jiang,  S.D.  Bader,  J.  Appl.  Phys.  93  (2004)  6507. 

[10]  W.  Hahn,  M.  Loewenhaupt.  Y.Y.  Huang,  G.P.  Felcher, 
S.S.P.  Parkin,  Phys.  Rev.  B  52  (1995)  16041. 

[11]  C.S.  Arnold,  D.P.  Pappas,  A.P.  Popov,  Phys.  Rev.  Lett.  83 
(1999)  3305. 

[12]  G.P.  Felcher,  J.  Appl.  Phys.  87  (2000)  5431. 

[13]  Ph.  Bauer,  M.  Sajieddine,  C.  Dufour,  K.  Cherifi, 
G.  Marchal,  Ph.  Mangin,  Europhys.  Lett.  16  (1991)  307. 

[14]  F.  Itoh,  N.  Nakamura,  H.  Sakurai,  J.  Magn.  Magn. 
Mater.  126  (1993)  361. 

[15]  N.  Ishimatsu,  H.  Hashizume,  S.  Hamada,  N.  Hosoito,  C.S. 
Nelson,  C.T.  Venkataraman,  G.  Srajer,  J.C.  Lang,  Phys. 
Rev.  B  60  (1999)  9596. 

[16]  H.  Sano,  H.  Hashizume,  H.  Okuda,  N.  Hosoito,  Jpn.  J. 
Appl.  Phys.  41  (2002)  103. 

[17]  N.  Hosoito,  H.  Hashizume,  N.  Ishimatsu,  I.-T.  Bae,  G. 
Srajer,  J.C.  Lang,  C.K.  Venkataraman,  C.S.  Nelson, 
Jpn.  J.  Appl.  Phys.  41  (2002)  131. 

[18]  N.  Hosoito,  H.  Hashizume,  N.  Ishimatsu,  J.  Phys. 
Condens.  Matter  14  (2002)  5289. 


